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ABSTRACT 


The  C12(d,n)N13  react  ion  has  been  studied  at  a  deuteron  bom¬ 
barding  energy  of  6  MeV0  Angular  distributions  of  the  neutrons  leading 
to  the  ground 8  2a36j  and  the  unresolved  3 » 51-3 # 56  MeV  states  were 
obtainedo 

The  ground  (J^  =  1/2”),  2„367  (J77  =  l/2+)  and  the  unseparated 
3o5 6  (J71  =  5/2+)  MeV  states  were  analyzed  with  the  plane  wave  Born 
approximation o  An  analysis  of  the  separated  3  <>51  (J77  =  3/2”)  and  3 <>56 
MeV  states  was  attempted0  The  ground  state  was  analyzed  using  the  dis¬ 
torted  wave  Born  approximation 

Spectroscopic  factors  were  extracted  for  the  four  levels  of  N130 
The  width  of  the  unseparated  doublet  in  N13  at  3o51-3<>56  MeV  was  found 
to  be  about  71  keV0  There  was  no  evidence  of  spin  dependent  structure 
for  this  nucleus  at  this  energy Q 

An  appendix  describes  the  conversion  of  two  elastic  search  codes 
along  with  an  outline  of  the  theory  relevant  to  (d,n)  reactions 0 
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I .  INTRODUCTION 


a.  Motivation  for  the  Experiment 

Nuclear  stripping  reactions  (G1  63)  have  provided  physicists 
with  an  extremely  useful  tool  for  studying  the  properties  of  nuclei „ 

This  is  due  to  the  fact  that  only  a  few  nucleons  are  actively  involved 
in  the  reaction 9  as  well  as  the  fact  that  hound  nuclear  states  can  be 
reached 0  From  physical  considerations  it  is  apparent  that  the  cross 
section  will  involve  the  overlap  of  wave  functions  describing  the  ini¬ 
tial  system  with  those  describing  the  final  system 0  Because  of  the 
few  degrees  of  freedom  involved  such  an  overlap  measures  the  degree 
to  which  the  passive  nucleons  occupy  the  same  configuration  in  the 
initial  and  final  states 0 

Furthermore  the  fact  that  there  are  only  three  active  particles 
involved  enables  an  unambiguous  determination  to  be  made  of  the  angu¬ 
lar  momentum  of  the  captured  particle 0  Conservation  of  angular  momen¬ 
tum  imposes  a  restriction  on  the  spin  state  of  the  final  nucleus „ 
Explicitly  we  must  satisfy 

Jj.  +  +  1/2  =  Jf,  »i»f  =  (-D*p  (1-1) 

for  a  (d,n)  reaction,,  is  the  angular  momentum  of  the  captured  proton, 

and  Jj_,  7T ^  and  ,  TTf  are  the  spins  and  parities  of  the  initial  and  final 
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nuclei  respectively . 

The  angular  distribution  of  the  emitted  neutron  is  frequently 
characteristic  of  the  value  of  and  if  Jj_  is  known  this  allows  close 
limits  to  be  set  on  J^..  At  best,  when  the  target  nucleus  has  zero 
spin,  the  final  spin  is  limited  to  J^.  =  l  ±  1/2.  The  remaining  ambi¬ 
guity  can  only  be  resolved  by  polarization  or  y-rey  correlation  exper¬ 
iments  0 

Lee  and  Schiffer  ( Le  6b)  have  recently  observed  for  (d,p) 
reactions  a  consistent  difference  in  the  structure  of  angular  distri¬ 
butions  of  states  populated  by  the  same  £-value  neutrons  but  differing 
in  Jfo  The  effect  was  most  evident  for  £n  =  1.  For  these  transitions 
and  deuteron  energies  between  7  and  12  MeV,  40  <  A  <  65 ,  and  spectro¬ 
scopic  factors  greater  than  about  0.1  (see  section  IV  d)  the  =  1/2 
states  exhibit  a  sharp  minimum  somewhere  in  the  range  90°  <  6  <  1^5° 
but  Jf=  3/2  states  do  not.  The  effect  decreases  with  decrease  in 
energy.  Experiments  with  light  nuclei  do  not  seem  to  show  this  trend. 

The  purpose  of  this  experiment  was  then  to  investigate  the 
presence  of  this  effect  in  light  nuclei,  in  particular  with  the 
C12(d,n)  reaction  to  the  ground  state,  J77  =  1/2  ,  and  to  the  second 
excited  state,  J77  =  3/2“.  It  was  necessary  to  split  the  doublet  in 
N13  at  3«51  -  3o56  MeV  excitation  in  order  to  extract  the  minor  com¬ 
ponent  (Figure  1-1)  „  A  bombarding  energy  of  6  MeV  provided  optimum 
spectrometer  resolution.  The  procedure  and  success  of  resolving  the 
doublet  is  described  in  section  III  b.  The  existence  of  an  excellent 
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fast  neutron  time-of -flight  spectrometer  made  this  work  at  all  possible* 

12,  .  13 

b»  Survey  of  the  Literature  on  the  C  (d.n)N  Reaction 

Beginning  in  191+0  and  up  to  1951  reaction  products  of  the  deu- 
teron  bombardment  of  carbon  had  been  studied  (Bo  1+9),  and  analysed 
assuming  that  the  compound  nucleus  was  the  dominant  source  of  these 
products*  Using  0.7  to  1,9  MeV  deuterons  and  a  proportional  counter, 
Bonner  (ibid,)  however  already  observed  that  the  ratio  of  proton  to 
neutron  yields  was  much  too  high  when  compared  with  Coulomb  penetra¬ 
bilities  0 

In  the  same  year  Grosskreutz  (Gr  1+9)  with  10  MeV  deuterons  frcm 

the  Washington  University  cyclotron  and  photoplates  first  reported  a 

1  3 

state  at  3.1+8  ±  0.12  MeV  excitation  in  J  .  Van  Patter  (Va  1+9)  subse¬ 
quently  associated  the  first  three  levels  of  N13  with  those  of  the 
mirror  nucleus  C13, 

In  a  series  of  experiments  at  Liverpool  in  1951  El-Bedewi, 
Middleton,  and  Tai  (El  51)  first  reported  the  angular  distribution  of 
separated  neutron  groups,  using  8  MeV  deuterons  and  photographic  emul¬ 
sions.  They  interpreted  the  neutrons  feeding  the  ground  state  (Q  = 
-0.21+  MeV),  2.38  ±  0.05  MeV  state  and  3.53  ±  0.05  MeV  state  of  N13  by 
Newn’s  stripping  theory,  assigning  spins  of  3/2"  or  1/2"  (  £p  =  1,  R  = 
5.5  f),  l/2+  (£p  =  0,  R  =  6.3  f),  and  5/2+  or  3/2+  Up  =  2,  R  =  6.3  f), 
respectively.  Butler  fits  (Mi  53)  yielded  radii  one  to  two  fermis 
smaller  (the  Gamow  radius).  However  for  a  radius  consistent  for  all 
the  states  Middleton  et  al.  found  that  Butler  theory  predicted  the 
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wrong  £~value  for  the  first  excited  state .  They  therefore  concluded 
that  for  this  particular  Q-value  and  their  deuteron  energy  Butler 
theory  is  incorrect;  prohahly  due  to  the  singular  behaviour  near  the 
binding  energy  of  the  deuteron .  They  found  other  discrepancies  canmon 
to  this  experiment o  These  were  an  experimental  forward  minimum  which  is 
not  predicted  for  £p  =  1  and  the  forward  maximum  for  the  3 » 53  MeV 
group  where  £  is  not  zero,  which  is  characteristic  of  high  negative 
Q  values o  If  one  allows  some  arbitrariness  in  the  radius  then  £  is 
not  uniquely  determined  in  this  latter  case.  Comparison  with  the 
mirror  nucleus  C13  and  the  shell  model  confirmed  the  ^i/2  an<^  ^1/2 
nature  of  the  ground  and  first  excited  states,  respectively.  Concern¬ 
ing  the  second  excited  state,  Jackson  and  Galonsky  (Ja  51)  first  re¬ 
ported  it  as  a  doublet  from  elastic  proton  scattering,  consisting  of 
a  ^3/2  s^a^e  a"k  3<>503  MeV  and  a  D5/2  state  at  3.550  MeV.  Rotblat 
(Ro  51)  confirmed  two  similar  states  in  C13  at  3.683  and  3o88U  MeV 
with  spins  or  P3/2  and  D3/2  or  D5/2»  respectively. 

Wilkinson  (Wi  55)  made  use  of  the  unique  property  of  N1 3  that 
its  excited  states  are  all  unbound  to  proton  emission  to  measure  total 
C12(d,n)N13  cross  sections  to  the  ground  state.  He  passed  a  20  MeV 
deuteron  beam  through  polyethylene  foils  stacked  between  sheets  of 
aluminum.  The  residual  0+  activity  plotted  against  time  and  extra¬ 
polated  back  to  when  the  beam  was  stopped  yielded  the  total  cross 
section.  He  found  a  smooth  excitation  function.  This  suggested  little 
compound  nucleus  formation,  although  at  the  higher  energies  the 
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aluminum  foils  represented  some  500  kilovolts  thickness.  The  observed 
total  cross  section  was  100  mb  at  =  8  MeV. 

Bonner  et  al.  (Bo  56)  studied  the  (d,n)  and  (d,p)  reactions 
from  C12  at  1.8  to  6.1  MeV  bombarding  energy  in  an  effort  to  extract 
the  direct,  compound,  and  interference  contribution  in  this  transition 
region.  The  deuteron  source  was  the  6  MeV  Rice  Institute  Van  de 
Graaff .  The  proton  detector  was  a  scintillation  counter  while  the 
neutron  detector  was  both  a  long  counter  and  an  energy-sensitive  de¬ 
tector  made  up  of  four  spherical  anthracene  scintillators.  With  a 
155  yg/cm2  target  the  cross  section  at  Ed  =  3.64  MeV  for  the  ground  state 
neutron  group  was  found  to  be  14  mb/sterad  at  30°.  Resonances  were 
found  at  deuteron  energies  of  2.47,  2.67,  2.99,  3.39,  4.00,  4.6.  4.8, 
5.34,  and  5.64  MeV.  Bonner  found,  by  comparison  with  Butler  theory 
and  by  comparing  total  with  first  excited  state  excitation  functions, 
that  the  neutron  (proton)  group  feeding  the  first  excited  state  of 
N1 3 ( C 1 3 )  do  so  more  by  stripping  than  those  feeding  the  ground  state. 
Bonner  points  out,  however,  that  this  may  be  unique  for  the  deuteron 
energy  used  ,  if  in  the  stripping  picture  one  assumes  that  the  disinte¬ 
gration  particle  has  about  the  same  velocity  as  the  incident  deuteron. 

By  writing  the  cross  section  as  a  coherent  sum  of  a  direct  (background) 
contribution  plus  a  resonant  compound  nucleus  contribution  and  sepa¬ 
rated  in  phase,  Bonner  found  that  for  the  (d,p)  reaction  at  4  MeV  and 
30°,  the  direct  contribution  was  63%.  Of  the  remaining  37%,  33%  was 
interference  and  h%  was  due  to  the  resonant  term  alone.  In  this  way 
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he  found  angular  distributions  symmetric  about  90°  and  extracted  spins 
of  Nll+#  consistent  with  other  work0 

Benenson  et  al 0  (Be  55)  made  a  comparison  of  the  differential 
cross  sections  for  the  reactions  C12(dsp)C12  and  C12(dgn)N13  at  bom¬ 
barding  energies  of  2o68  MeV  and  3o26  MeV0  Using  solid  and  gas  targets 
and  a  gas  recoil  neutron  spectrometer  they  obtained  absolute  cross 
sections  and  fitted  these  to  Butler  theory 0  For  the  (d9p)  reaction 
they  first  subtracted  an  isotropic  background 0  In  view  of  the  large 
possible  interference  this  was  not  a  good  idea0  In  contrast  with  the 
findings  of  Middleton  et  al Q  (Mi  53 ) s  the  resulting  agreement  with 
Butler  theory  for  the  ground  state  was  quite  good0  They  then  compared 
the  ratios  of  the  (d9n)  cross  sections  with  the  ratios  of  the  respec¬ 
tive  kinetic  factors  from  plane  wave  stripping  theory 0  They  concluded 
that  the  reasonable  agreement  was  justified 0  Of  the  effects  neglected,, 
effects  involving  the  deuteron  appear  in  both  reactions 0  The  outgoing 
proton  had  an  energy  over  5  MeVs  well  above  the  Coulomb  barrier 0 
Furthermore  the  nuclear  scattering  phase  shifts  of  the  liberated  par¬ 
ticles  s  appreciable  only  for  angular  momenta  less  than  two8were  found 
to  be  similar o 

Calvert  et  al0  (Ca  57)  bombarded  carbon  with  9  MeV  deuterons 
and  found  agreement  with  Butler  theory  consistent  with  the  work  of 
Middleton  et  al0  (Mi  53)  *  but  with  a  radius  parameter  of  (l<>7  + 

1 0 22A1 / 3 )  fermis2  which  is  the  Gamow  radius »  In  addition  they  meas¬ 
ured  the  absolute  cross  sections  for  the  three  neutron  groups ,  includ- 
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ing  the  unresolved  doublet 0 

McEllistrem  (Me  58)  made  an  extensive  analysis  of  the  C12  +  d 
reactions  showing  that  the  resonances  in  N14*  were  of  compound  nucleus 
origin0  He  found  corrected  reduced  widths  for  the  (d,p)  reaction  by 
canparing  only  the  higher  angular  momentum  components  of  an  extracted 
reaction  amplitude  with  that  of  the  Butler  amplitude „  He  concluded 
that  the  data  could  best  be  represented  by  a  stripping  amplitude  plus 
more  than  one  resonant  compound  nuclear  amplitudes 0 

Elwyn,  Kane,  0fer8  and  Wilkinson  (El  59)  »  did  the  first  time- 
of “flight  experiment  with  C12(d,n)N13,  at  lU  energies  frcm  E^  =  I0U5 
MeV  to  Ed  =  2o95  MeV0  They  found  that  the  direct  reaction  was  domi¬ 
nant  at  all  the  bombarding  energies 0  Fluctuation  in  the  best  fit 
Butler  radii  were  most  likely  due  to  the  large  direct-compound  nuclear 
interference  term0  They  attempted  to  remove  this  term  by  averaging 
the  cross  section  at  the  first  peak  over  several  deuteron  energies 0 
In  this  way  the  effects  of  the  excitation  energy  (equivalently  the 
bombarding  energy)  sensitive  compound  state  would  be  eliminated,.  The 
resulting  fit  to  Butler  theory  was  quite  good  in  the  forward  hemi¬ 
sphere  0  That  the  compound  nuclear  contribution  itself  is  small  was 
evidenced  in  the  smoothness  of  the  total  integrated  cross  section.  It 
has  been  pointed  out  (Wi  58),  in  fact,  that  in  low  Q  stripping  with 
low  energy  deuterons  direct  effects  should  be  large  on  simple  physical 
grounds  (see  also  Gi  6l).  Consideration  of  the  binding  energy  of  the 
proton  in  N13,  l„9Ul  MeV,  and  the  neutron  in  C12,  18,722  MeV  (Ho  6l), 
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led  to  the  conclusion  that  projectile  stripping  was  the  main  mechanism,, 
Therefore  Coulomb  and  nuclear  effects  have  to  be  included  to  obtain 
overall  agreement ® 

James  (Ja  6l)  has  investigated  the  C12(d,n)N13  reaction  below  2 
MeV  bombarding  energy 0  He  found  agreement  consistent  with  compound 
nucleus  formation  but  concluded  that  charge  symmetry  in  the  (d,p) 

-  (d,n)  reactions  was  not  conserved  because  of  the  presence  of  an  a4 
coefficient  in  the  (d,n)  case® 

Kolonko  et  ale  (Ko  63)  bombarded  C12  with  13  MeV  deuterons  and 
analysed  the  neutrons  leading  to  the  ground s  first  and  the  unresolved 
second  and  third  states  using  Butler  theory ®  Their  inconsistent 
radii  were  attributed  to  the  singularity  in  the  theory,  i»ec  the  state 
of  the  final  nucleus  is  an  energy  state  close  to  the  separation  energy 
of  the  last  nucleon  (Bh52)®  They  also  found  that  Coulomb  corrections 
to  the  Butler  formula  produced  a  stronger  variation  of  the  reduced 
width  with  rQ  than  without  this  correction® 

Jeronymo  et  al®  (Je  63)  performed  the  C12(d,n)C13  and 
C12(d,d)C12  reactions  at  bombarding  energies  from  0®9  to  5  MeV „  They 
found  that  there  is  a  large  resonant  contribution  in  deuteron  reactions 
with  C12  even  at  5  MeV®  They  point  out  that  since  in  this  nucleus  the 
p  subshell  is  filled,  it  is  possible  that  due  to  lack  of  loosely 
coupled  surface  nucleons  the  incident  particle  tends  to  form  a  compound 
system®  Hence  the  relative  contribution  of  resonance  mechanism  compared 
to  direct  interactions  gets  enhanced® 
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Hodgson  (Ho  63)  analysed  the  (d,n)  data  of  Elwyn  et  al ,  (El  59) 
with  a  zero  range  distorted  wave  theory  and  an  arbitrary  cut-off 
radius o  He  found  several  sets  of  optical  potential  parameters  that 
accounted  for  the  overall  distribution,.  The  absolute  cross  sections 
however,  as  in  this  experiment,  were  found  particularly  sensitive  to 
the  choice  of  potentials 0  Ohlsen  and  Shamu  (01  63)  have  recently 
scattered  6  MeV  deuterons  on  carbon0  Optical  potential  parameters 
from  their  data  would  no  doubt  prove  useful,  though  some  of  the  param¬ 
eters,  in  particular  the  absorbing  part  of  the  potential,  would  be  more 
accurately  fixed  by  reaction  data0  The  (d,d)  excitation  function  at 
6  MeV  lies  on  the  edge  of  the  prominent  15=09  MeV  level  and  on  the  less 
prominent  15=^4  MeV  levels  of  N14*,  This  would  further  affect  a  direct 
comparison  of  elastic  with  reaction  parameters, 

Fulbright ,  Robbins,  West,  and  Saylor  (Fu  66)  bombarded  carbon- 
12  at  energies  near  the  U  MeV  resonance  in  They  found  that 

at  small  angles  the  ciinimum  value  of  the  resonant  cross  section  to  the 
ground  state  of  N1 3  consistent  with  their  analysis  is  5  mb/sr0  The 
maximum  resonant  contribution  for  the  (d,p)  case,  on  the  other  hand, 
has  been  shown  by  Bonner  et_  al „  (Bo  56)  to  be  no  greater  than  1  mb/sr, 
While  Bonner  shows  that  in  the  (d,p)  case  to  the  ground  state  of  C13 
the  resonant  term  is  small  compared  to  the  resonant— direct  interference 
term,  Fulbright  et  al,  find  that  there  is  little  interference  in  both 
the  (d,n)  and  (d,p)  reactions  to  the  first  excited  states  of  N13  and 
C13,  respectively,  at  the  same  bombarding  energy. 
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The  foregoing  survey  has  shown  that  a  proper  treatment  of  the 
C12(d,n)N13  reactions  requires  a  contribution  from  compound  nucleus 
formation.,  This  is  probably  due  to  the  closed  shell  structure  of  C12 
and  the  resulting  lack  of  available  single  particle  states «,  Conclu¬ 
sions  drawn  from  the  present  experiment  could  therefore  be  invalidated 
by  such  unfortuitous  interference  as  studied  in  the  above  section, 
although  (d,p)  experiments  to  the  ground  state  of  C13  at  higher  ener¬ 
gies  indicate  that  this  is  unlikely „ 


1  1 
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II.  THE  EXPERIMENT 

a.  Experimental  Arrangement 

The  C12(d,n)N13  reaction  was  performed  with  the  University  of 
Alberta  6  MeV  model  CN  *)  Van  de  Graaff  accelerator,  utilizing  the 
time-of-f light  method  with  top  terminal  pulsing  and  a  75"  Mobley  com¬ 
pression  system  (Da  66).  The  beam  transport  system  is  illustrated  in 
Figure  2-1. 

Gas  is  ionized  in  the  source  bottle  by  a  125  me  source  oscil¬ 
lator  o  The  stripped  electrons  are  attracted  to  the  positive  probe 
while  the  positive  ions  are  attracted  to  an  opening  called  the  canal . 
The  helical  motion  of  the  electrons  is  due  to  the  source  magnet. 

This  also  enhances  ionization.  Since  the  net  charge  in  the  resulting 
plasma  is  zero,  the  potential  between  the  probe  and  canal  is  confined 
to  the  Crookes  dark  space  near  the  canal. 

The  terminal  analysis  magnet  acts  as  a  mass  separator.  Focus¬ 
ing  is  accomplished  by  a  three  element  electrostatic  Einzel  lens,  the 
self  focusing  column,  and  a  quadrupole  doublet. 

Beam  pulsing  is  achieved  by  two  sets  of  parallel  deflection 
plates  perpendicular  to  each  other.  These  are  driven  out  of  phase  by 

1 )  High  Voltage  Engineering  Corporation,  Burlington,  Massachusetts, 
U.SoAo 
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FIGURE  2-1  BLOCK  DIAGRAM  OF  BEAM  TRANSPORT  SYSTEM 
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a  1  me  RoFo  oscillator „  The  resulting  elliptical  motion  of  the  team 
over  an  aperture  produces  the  desired  pulsing 0 

Energy  (mcmentum)  analysis  of  the  beam  is  done  by  the  analyzing 
magneto  The  switching  magnet  serves  to  deflect  the  beam  to  the  appro¬ 
priate  experimental  port» 

The  10  me  RcF0  deflection  system  2)  deflects  the  individual 
beam  bursts  for  compression  in  the  Mobley  magneto  Compression  is  real¬ 
ized  by  virtue  of  the  same  time  required  for  the  two  ends  of  the  burst 
to  traverse  a  different  distance »  The  original  10  ns  bursts  are  com¬ 
pressed  in  this  manner  to  0„5  ns0  The  resulting  average  beam  current 
is  of  the  order  of  2  y  amps0 
bo  The  Neutron  Time-of-Flight  Spectrometer 

The  time-of-f light  spectrometer  is  illustrated  in  Figure  2-2 « 

The  stop  signal  is  derived  frem  a  3  cm  long  cylindrical  beam  sensing 
capacitor  located  30  cm  before  the  target 0  The  start  signal  ccraes  from 
the  neutron  detector 0  This  is  a  cylindrical  Ne  213  3)  3oU5"  x  0<,75" 
liquid  scintillator  enclosed  in  quartz ,  It  is  coupled  to  a  Phillips 
XPIO^O  photomultiplier o  A  Naton  phosphor  4)  coupled  to  an  RCA  8575 
photomultiplier  serves  as  the  monitor 0 

The  main  detector  was  mounted  on  a  remote  controlled  cart  at  the 
end  of  a  six  meter  boom0  Neutrons  from  carbon  in  the  system  were 

2)  See  (Da  66)  and  references  therein  „ 

3)  Nuclear  Enterprises  Ltd.,  Winnipeg,  Manitoba,, 

4)  Nash  and  Thompson  Ltd„  ,  Hookrise  South,  Tolworth,  Surrey,  England,, 
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FIGURE  2-2  NEUTRON  TIME  OF  FLIGHT  SYSTEM 
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shielded  by  concrete,  water  and  paraffin  up  to  within  30  inches  of  the 
target  chamber  and  thereon  by  two  inches  of  lead. 

The  time-to- amplitude  converter  was  based  on  a  design  by  D#L# 
Wieber  (Wi  63)  using  capacitance  charging#  The  time  base  was  about 
300  nso  The  output  of  the  TAC  was  fed  to  a  TMC  102U  channel  pulse 
height  analyzer  (ADCA)0  The  linear  channel  from  dynode  11  of  the 
XPlOl+O  accepts  only  neutrons  above  a  threshold  energy,,  This  channel 
is  also  used  for  pulse  "walk"  correction#  Walk  is  broadening  of  the 
peaks  in  the  spectrum  due  to  the  difference  in  observed  timing  caused 
by  the  different  recoil  proton  pulse  heights#  This  channel  carries  the 
pulse  height  information#  Figure  2-2a  shows  an  uncorrected  3-D  spectrum# 
Cross  over  timing  is  used  in  the  monitor  start  channel  as  this 
is  less  susceptible  to  walk#  Leading  edge  timing  however  gives  better 
inherent  resolution#  Since  the  neutron  channel  is  corrected  for  walk, 
leading  edge  timing  is  therefore  used#  The  stop  signal  provides  a 
pulse  with  the  centroid  of  the  passing  beam  burst  defined  by  the  cross¬ 
over  as  the  induced  charge  swings  from  positive  to  negative#  Cross-over 
timing  is  therefore  used  in  the  stop  channel# 

An  acceptance  pulse  from  the  gamma  ray  eliminator  (Re  66)  in 
coincidence  with  a  prompt  signal  from  the  linear  channel  enters  the 
coincidence  interrupt  logic  unit  (Ge  66)#  Here  the  coincidence  pulse 
has  to  pass  the  two  ADS  Busy  inhibit  gates  and  then 

l)  energize  the  two  ADC  units  via  the  coincidence  interrupt 
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2)  in  triple  coincidence  with  the  two  ADC  Address  Reset  signals 
energize  ADC  transmission  to  the  computer. 

Walk  correction  is  accomplished  by  dividing  the  linear  signal 
into  32  amplitude  bins  and  finding  the  centroid  of  each  bin.  The  dis¬ 
tance  in  channels  from  a  channel  in  the  reference  bin  is  the  shift. 

The  resulting  shift  table  is  read  back  in  and  the  correction  applied 
when  in  the  NTOF  mode. 

The  overall  resolution  was  about  0.9  ns  as  compared  with  a  pre¬ 
vious  figure  of  1.25  ns  (Da  66)  using  phototube  saturation  to  correct 
for  walk.  The  major  improvement  in  the  present  method  is  at  the  base 
of  the  peaks.  Using  AE/E  =  2At/t  we  have  a  resolution  of  about  45  keV 
for  5  MeV  neutrons  traversing  a  6  meter  flight  path. 
c.  Efficiency  of  the  Neutron  Detector 

The  relative  efficiency  of  the  neutron  detector  was  determined 
by  a  separate  experiment  using  the  T(p,n)He3  reaction.  The  energy 
calibration  of  the  accelerator  was  previously  done  by  W.G.  Davies 
(Da  66)  using  the  Li7(p,n)Be7  reaction.  The  neutrons  were  detected  by 
the  main  neutron  detector  and  by  a  McKibben  long  counter  with  paraffin 
moderator  (AJ.  60).  The  BF3  electronics  are  illustrated  in  Figure  2-3. 
The  integrated  count  rate  was  corrected  for  counting  efficiency  and 
used  as  the  normalization  for  the  relative  efficiency  curve.  The 
absolute  efficiency  was  obtained  from  an  associated  particle  time-of- 
f light  experiment  (Bu  66)  using  the  D(d,n)He3  reaction  in  conjunction 
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BF3  ELECTRONICS  FOR  RELATIVE  EFFICIENCY  MEASUREMENT 
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The  agreement  of  the  efficiency  versus  neutron  energy  with  the 
theoretical  efficiency 


(2-1) 


is  shown  in  Figure  2-b0  In  this  expression 
En  is  the  neutron  energy 0 

nj|  is  the  number  of  hydrogen  atoms  per  c0c0  calculated  for 
Ne  213  as  U,888  x  1022  hydrogen  atoms/c ,  c,  (Nu  65, 

Re  66) 0 

on_p  is  the  total  neutron-proton  cross  section  (Ba  57 ) « 
x  is  the  detector  length,  in  this  case  1,905  cm, 

Agreement  with  theory  is  not  very  good.  As  a  result  a  cubic  inter¬ 
polation  was  performed  between  smoothed  experimental  points.  This 
was  fitted  asymptotically  to  the  n-p  cross  section.  There  seems  to 
be  no  obvious  reason  for  this  lack  of  agreement, 

A  further  small  correction  was  applied  to  the  forward  angle 
ground  state  neutron  group  due  to  an  upper  cut-off  by  the  single 
channel  analyzer  at  a  neutron  energy  of  5«35  MeV ,  This  was  done  with 
a  two-dimensional  Be9(d,n)B10  spectrum  which  presents  pulse  height  vs, 
neutron  energy  vs,  recoil  proton  energy.  Extrapolation  to  the  cut-off 
energy  was  made  using  the  expression  for  pulse  height  vs,  proton  energy 
from  Gove  et  al ,  (Go  6l) 


Lp  =  0.215  Ep  +  0.028  Ep2  (0  <  E  <  8  MeV) 


(2-2) 
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FIG-  2-4  EFFICIENCY  OF  THE  NEUTRON  DETECTOR 
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for  liquid  NE  213 0 

As  a  consistency  check  for  future  efficiency  measurements  a  thick 
yield  was  taken  from  the  Mn^^(dsn)Fe^^  reaction,,  A  yield  with  the  gamma 
eliminator  in  is  shown  in  Figure  2-5 „ 

Figure  2-6  shows  the  yield  as  a  function  of  neutron  energy  for 
both  the  low  and  high  ranges 0  At  the  time  of  the  experiment  the  dynamic 
range  of  pulse  heights  accepted  by  the  side  channels  had  to  be  divided 
into  two  regions,  corresponding  to  neutron  energies  approximately  0o5 
to  5  MeV  and  1  to  10  MeV0 
d0  Targets 

Carbon  targets  were  initially  prepared  by  arc  evaporation  as 
described  by  G0  Dearnaley  (De  6l)„  However  it  was  difficult  to  obtain 
a  thickness  over  20  yg/cm2,, 

Thin  carbon  targets  were  then  prepared  by  smoking  clean  0,005" 
tantalum  blanks  in  a  benzene  flame,,  Figure  2-7  shows  the  benzene 
burner  0 

The  thickness  of  the  target  was  measured  by  weighing  and  deter¬ 
mined  to  be  approximately  50  ±  5  yg/cm0  The  target  was  reweighed  after 
the  experiment  and  then  one-half  and  one-quarter  of  the  target  were 
weighed.  These  measurements  were  consistent  with  the  previous  deter¬ 
mination  of  50  yg/cm20 
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III-  DATA  analysis-aquisition  and  reduction  of  data 

a*  Procedure 

Angular  distributions  of  the  neutrons  from  the  C12(d,n)N13 
reaction  were  taken  at  twenty  angles;  in  5°  steps  from  0°  to  65° 
and  in  10°  steps  from  75°  to  IU50,  except  for  two  angles  at  35°  and 
55 °°  The  detector  was  placed  at  a  distance  of  six  meters  from  the 
target  and  the  monitor  was  located  at  1„5  meters  and  -30° .  The  live 
times  of  the  monitor  channel  and  the  main  detector  channel  were  re¬ 
corded  and  these  were  used  to  correct  the  yield,  The  analyzing  magnet 
NMR  averaged  out  to  29*296  ±  0,001  megacycles  which  corresponds  to  a 
deuteron  bombarding  energy  of  5*99^  MeV, 

A  typical  spectrum  is  shown  at  0°  in  Figure  3-1*  The  raw  data 
was  fitted  with  a  nonlinear  least  squares  Gaussian  fitting  program 
(Da  66)  which  subtracted  background  of  a  given  shape ,  A  linear  back¬ 
ground  was  chosen  as  this  least  affected  the  area  of  the  Gaussian* 

The  Gaussian  in  most  cases  underestimated  the  area  by  about  3%  as  there 
was  no  skewness  compensation  and  the  actual  peaks  were  more  pointed, 
although  a  variable  standard  deviation  compensated  somewhat  for  this* 
The  fitted  areas  were  nevertheless  used  as  the  computer  was  more  con¬ 
sistent*  Furthermore  the  monitor  spectrum  was  also  fitted  and  this 
tended  to  cancel  an  underestimate  of  the  peak  area. 
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Absolute  cross  sections  were  calculated  from  the  expression 

C  =  o(6)  •  dfi  •  N  •  /  Idt  •  e  (3-1 ) 

where 

C  is  the  total  number  of  counts  under  the  peak0 
o(0)  is  the  cross  section  in  millibars /steradian. 
dft  is  the  solid  angle,  which  in  this  case  was  1 0666  *  10 
steradians  for  a  detector  distance  of  6  meters „ 

N  is  the  number  of  target  nuclei/cm2  calculated  to  be 
(2o59  ±  0o2 6)  x  10 18  target  nuclei/cm2 0 
j  Idt  is  the  integrated  beam  current  on  the  targets 
e  is  the  absolute  efficiency  of  the  neutron  detector  calcu¬ 
lated  as  in  section  II  b  „ 

The  beam  current  was  integrated  with  an  Elcor  Model  A309A  current 
integrator,,  The  monitor  was  used  for  normalization  by  finding  the 
average  ratio  of  monitor  counts  to  integrated  beam  current  for  several 
runs  where  double  pulsing  was  minimal,,  Double  pulsing  is  caused  by 
the  elliptically  swept  beam  at  the  top  terminal  passing  over  the  aper¬ 
ture  twice  in  each  cycle „  This  ratio  then  gave  the  corrected  inte¬ 
grated  current  for  all  the  runs „ 
b0  The  3 o 51-3 0^6  MeV  Doublet  in  N1  3 

The  two  states  involved  were  essentially  unresolved,  although 
the  presence  of  the  weakly  excited  3° 51  MeV  state  was  clearly  notice¬ 
able  o  The  nonlinear  Gaussian  least  squares  program  was  used  in  an 
effort  to  split  this  doublet „  No  previous  neutron  work  has  succeeded 
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in  splitting  the  two  states,  and  this  is  further  complicated  by  the 
fact  that  while  these  two  states  are  only  50  keV  apart ,  they  are  respec¬ 
tively  63  and  6l  keV  wide  (a  j  59 )  c 

There  were  two  methods  of  approaching  the  problem .  One  was 
to  obtain  the  peak  means  from  a  time  calibration  line  extracted  from  a 
Beryllium  reference  spectrum.  The  Beryllium  peaks,  however,  were  situ¬ 
ated  at  the  high  energy  end  of  the  kicksorter  and  the  known  nonlinear¬ 
ity  of  the  ADC  would  render  extrapolation  meaningless. 

Alternatively  one  could  assume  that  the  3.51  MeV  state  is  weakly 
excited  with  respect  to  the  3.56  MeV  state  and,  assuming  two  Gaussians, 
search  for  the  best  peak  means.  The  weak  excitation  of  the  3*51  MeV 
level  is  evidenced  in  the  mirror  levels  in  C13  at  3*68  and  3.85  MeV 
where  McGruer  et  al.  (Me  55)  found  that  a 3 a 85 5 (  © ) /cr 3  535(0)  is  9  =  5. 

From  elastic  proton  scattering  on  C12  Jackson  and  Galonsky  (Ja  51) 
found  a  ratio  of  6.3  for  the  reduced  widths  of  the  3.56  to  3.51  MeV 
levels . 

With  this  in  mind  the  best  fit  mean  to  the  3.56  MeV  state,  keep¬ 
ing  the  standard  deviations  fixed,  was  taken  to  be  the  true  mean  of 
that  level.  The  position  of  the  3.51  MeV  level  was  fixed  by  finding 
the  number  of  channels  corresponding  to  the  known  kinematic  time  differ¬ 
ence  (Gr  65)  and  interpolating  for  the  curve  in  Figure  3-2  corresponding 
to  an  ADCA  baseline  setting  of  0.1.  This  technique  yielded  means  for 

f. 

the  3.51  MeV  state  consistent  with  a  visual  estimate.  The  integral 
linearity  of  the  analog  to  digital  converter  is  shown  in  Figure  3-3. 
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Since  these  two  states  have  almost  the  same  width „  this  parameter 
was  fixed  to  the  same  constant  in  the  fitting  procedure 0  The  shape  of  the 
extracted  angular  distribution  for  the  3o51  MeV  state  was  unaffected  by 
the  value  of  the  standard  deviation  chosen 0  This  indicates  that  what  was 
actually  extracted  was  probably  nothing  more  than  a  fairly  constant  por¬ 
tion  of  the  total  peak0  The  instrumental  width 9  including  target  thick¬ 
ness  s  was  much  less  than  the  total  observed  width 0  Therefore  further  im¬ 
provement  in  time  resolution  is  not  likely  to  yield  better  separation 0 
Application  of  a  more  realistic  fitting  program  (Te  66)  might  prove  more 
fruitful o  A  time  calibration  line  to  fix  the  peak  means  was  tried  but 
the  extrapolated  means  were  too  uncertain  to  be  of  any  use0 

It  was  not  possible  to  deduce  directly  the  width  of  the  2036Y  MeV 
state  as  there  were  no  nearby  peaks  from  which  to  extract  the  instru¬ 
mental  width o  With  neutrons  of  the  same  energy  from  the  Li7(p,n)Be7 
reaction  for  this  purpose,  a  width  of  ^9  keV  was  found0  This  is  not  in 
good  agreement  with  the  previously  reported  value  ( AJ  59)  of  32  keV 0 
The  combined  width  of  the  3 0 51*3 o 56  MeV  doublet  could  be  estimated  since 
here  the  instrumental  width  is  negligible  compared  to  the  widths  of  the 
two  components o  The  value  arrived  at  after  unfolding  an  instrumental 
width  extracted  from  the  N13  ground  state  neutron  group  was  71  keV0 
This  is  in  agreement  with  a  6l  keV  level  at  3o56  MeV  folded  in  with  a 
63  keV  level  at  3o51  MeV  and  one-sixth  as  intense 0 
c  o  Absolute  Cross  Sections 


Absolute  cross  sections  were  obtained  from  the  raw  data  using  a 
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computer  program  written  by  Mr0  T„B.  Grandy  which  corrected  the  meas¬ 
ured  intensities  for  neutron  attenuation  in  the  gold  backing  and  made 
the  necessary  laboratory  to  center  of  mass  trans format ion ,  The  overall 
normalization  was  obtained  using  the  procedure  of  section  III  a,  The 
efficiency  subroutine  was  rewritten  to  accomodate  a  cubic  interpola¬ 
tion  asymptotic  to  the  o  cross  section e  Tables  III-l  to  III-5  list 

n-p 

the  absolute  cross  sections  found  in  this  experiment. 

The  relative  errors  shown  in  the  tables  are  in  absolute  units. 
These  errors  are  from  the  Gaussian  fitting  program  and  from  the  effic¬ 
iency,  Not  shown  is  the  overall  normalization  error  arising  from 
target  thickness  and  beam  integration.  This  adds  another  10$  error 
giving  a  total  error  of  about  20$ . 

Table  III-3  contains  the  combined  doublet  in  N130  This  is  taken 
to  be  the  cross  section  for  the  3o56  MeV  component  for  the  reasons  out¬ 
lined  in  section  III  b0  Tables  III-H  and  III-5  contain  the  cross  sec¬ 
tions  for  each  component  of  the  doublet,  as  separated  by  the  Gaussian 
fitting  programo  The  total  errors  in  these  values  are  much  larger, 
that  for  the  3,51  MeV  level  is  probably  between  50$  and  100$, 
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ABSOLUTE  CROSS  SECTIONS  FOR  C12(d,n)N13 
TO  THE  GROUND  STATE  AT  ED  =  5.991*  MEV 


Angle 

Angle 

Cross  Section 

Relative 

(Lab) 

(C.  of  M.) 

(mb/sr ) 

Error 

0.0 

0.0 

20.96 

1.95 

5.0 

5.6 

21.35 

2.20 

10  ,0 

11.2 

23.06 

2.15 

15.0 

16.7 

19.1+2 

1.96 

o 

o 

o 

CVJ 

22.3 

18.30 

1.66 

25  o0 

27.9 

11.95 

l.lU 

30.0 

33.  k 

7.02 

0.70 

1+0.0 

M+.3 

2.86 

0.25 

1+5.0 

1+9.8 

1.76 

0.18 

50.0 

55.1 

2.11 

0.18 

60.0 

65.9 

3.2U 

0.27 

65 .0 

71.1 

3.99 

0.33 

75 .0 

81.5 

3.98 

0.30 

85  oO 

91.7 

2.86 

0.23 

95»0 

101.7 

2.71 

0.22 

105.0 

111.5 

2. 3b 

0.20 

115.0 

121.1 

2.88 

0.23 

125.0 

130.5 

3 .66 

0.30 

135.0 

139.8 

3.21 

0.35 

ll+5-0 

1U8.9 

3.37 

0.36 

Table  III  -  1 
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ABSOLUTE  CROSS  SECTIONS  FOR  C12(d,n)N13 
TO  THE  2o367  MEV  STATE  AT  Ep  =  5.99*+  MEV 
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I 

i 

t 

I 
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; 


Angle 

(Lab) 

Angle 
(C.  of  M. ) 

Cross  Section 
(mb/sr ) 

Relative 

Error 

0.0 

0.0 

57.37 

1+.50 

5.0 

5.8 

56.27 

5.09 

10  .0 

11.6 

1+1+.10 

3.52 

15.0 

17. 4 

25.71 

2.33 

20.0 

23.2 

18.01 

1.1+1 

25  oO 

29.0 

8.26 

0.72 

30.0 

3^.7 

5 .61+ 

0.1+9 

1+0.0 

1+6.1 

5.3U 

0.38 

1+5.0 

51.7 

5.01+ 

0.1+0 

50,0 

57.2 

5.33 

0.36 

60.0 

68.2 

1+.5U 

0.32 

65  oO 

73.6 

3.58 

0.25 

75  «0 

Bh.i 

3.1+6 

0,23 

85  oO 

9H.1+ 

3.20 

0.22 

95.0 

101+.1+ 

3.66 

0.26 

105.0 

llU.l 

3.83 

0.28 

115.0 

123  .6 

3.97 

0,28 

125.0 

132.7 

1+.00 

0.28 

135.0 

lUl.7 

3.13 

0.23 

1U5.O 

150.1+ 

2.81 

0.21 

Table  III  -  2 
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ABSOLUTE  CROSS  SECTIONS  FOR  C12(d,n)N13  TO 
THE  3.5 6  NEV  STATE  (UNSEPARATED)  AT  ED  =  5.994  MEV 


Angle 

(Lab) 

Angle 
(Co  of  M.) 

Cross  Section 
(mb/sr ) 

Relative 

Error 

OoO 

OoO 

1+1,65 

3.15 

5.0 

60I 

U5.87 

4.00 

10o0 

12  0  3 

1+0.53 

3  oil 

15  oO 

18.1+ 

30.52 

2.66 

20  .0 

24.5 

30.98 

2.35 

25  oO 

30  o5 

21.19 

1.76 

30.0 

36  0  5 

19.27 

1. 57 

40.0 

U8.U 

14.07 

0.95 

1+5  oO 

54.  3 

8.68 

0  067 

50  oO 

6O0I 

7.51 

0.49 

60o0 

71.4 

4.68 

0.33 

65  .0 

76.9 

3.35 

0.24 

75  <>0 

87.7 

3.55 

0.24 

85  oO 

98ol 

4.42 

0.31 

95o0 

IO80I 

6.28 

0.43 

105  oO 

117  0  7 

7.36 

0.53 

115  oO 

1260  9 

10.15 

0.76 

125 .0 

135.8 

7.10 

0.57 

135  "0 

144.3 

6.97 

0.6l 

145.0 

152  0  5 

6.15 

0.57 

Table  III  -  3 


ABSOLUTE  CROSS  SECTIONS  FOR  C12(d,n)N13  TO 
THE  3o56  MEV  STATE  (SEPARATED)  AT  ED  =  5.994  MEV 


Angle 

Angle 

(Lab) 

(C.  of  M. ) 

0.0 

0.0 

5.0 

6.1 

10.0 

12.3 

15  oO 

18.4 

20.0 

24.5 

25  oO 

30.5 

30.0 

36.5 

4o,o 

1*8.1* 

45.0 

54.3 

50.0 

60.1 

60.0 

71.4 

65  oO 

76.9 

75  oO 

87.7 

85  oO 

98.1 

95  «0 

108.0 

105  oO 

H7.7 

115.0 

126.9 

125.0 

135.8 

135.0 

144.3 

1U5.O 

152.5 

Cross  Section 
(mb/sr ) 

Relati' 

Error 

28.58 

2.23 

35.92 

3.18 

30.84 

2.42 

24.32 

2.17 

25.03 

1.93 

15.39 

1.32 

14.40 

1.20 

11.11 

0.77 

6.48 

0.52 

5.25 

0.35 

3.53 

0. 25 

2.49 

0.18 

2.98 

0.95 

3.24 

0.23 

4.18 

0.29 

6  0O6 

0.44 

6.29 

0.47 

5.39 

0.42 

i*  .97 

0.38 

4.86 

0.39 

Table  III  -  4 


ABSOLUTE  CROSS  SECTIONS  FOR  C12(d,n)N13  TO 
THE  3.51  MEV  STATE  (SEPARATED)  AT  Ep  =  3«51  MEV 


Angle 

(Lab) 

Angle 
(C,  of  M, ) 

Cross  Section 
(mb/sr ) 

Relative 

Error 

OoO 

0,0 

9  .i+5 

0.89 

5  .0 

6,1 

7.09 

0.83 

10  oO 

12,2 

6.71 

0.73 

15  oO 

18,3 

6,01 

0,72 

20,0 

24,  4 

4,66 

0,51 

25  ,0 

30.4 

4 ,16 

0,46 

30,0 

36,4 

3.62 

0.62 

4o.o 

48.3 

2.15 

0,23 

45  oO 

54,1 

1.97 

0.23 

50,0 

59.9 

1.59 

0.15 

6o,o 

71o2 

0.93 

0,10 

65.0 

76.7 

0,60 

0,09 

75.0 

87.5 

0.51 

0,09 

85  «0 

97o9 

0,92 

0,12 

95cO 

107  0  9 

1.51 

0,15 

105o0 

117.5 

0.92 

0,20 

115.0 

126.7 

1.76 

0,22 

125  oO 

135.6 

0.98 

0,21 

135  oO 

144.1 

0.98 

0.19 

145.0 

152.4 

0.95 

o,  16 

Table  III  -  5 
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IV u  DATA  ANALYSIS-INTERPRETATION  OF  DATA 


a0  Plane  Wave  Analysis 

The  angular  distribution  extracted  from  the  observed  levels 
of  N13  were  analyzed  with  the  plane  wave  Born  approximation  using  the 
programmed  Lubitz  (Lu  57)  calculation,,  In  the  case  of  capture  into 
unbound  levels  of  the  final  nucleus „  however9  it  is  very  difficult  to 
determine  the  Jl-value  of  the  captured  particle  (Ma  6o)„ 

This  difficulty  is  born  out  in  the  present  experiment  (Fig¬ 
ure  l+«l)o  For  the  unseparated  neutron  group  at  3o51"3o56  MeV  excita¬ 
tion  in  N13  L  =  1  and  Z  =2  fitted  equally  wello  If  a  radius 
P  P 

parameter  between  r^  =  3„5  f  and  r^  =  8,0  f  is  acceptable 9  consistent 
with  the  literature  (Ma  60 8  Da  52)  8  then  Z  =  l(r  =  3<>0  f)  is  excluded 

sr  i r 

in  favor  of  Jl  =  2  (r  =  5  „5  f )  „ 

Jr  J tr 

Plane  wave  fits  to  the  separated  components  of  the  3o51-3»56 
MeV  doublet  (Figure  h -2)  are  more  ambiguous  in  the  choice  of  £-value 0 
This  possibly  reflects  the  uncertainty  in  the  separation  procedure 
(Chapter  III)  rather  than  ambiguity  in  the  theory 0  No  conclusion  can 
be  drawn  on  the  best  H-value0  The  same  is  true  of  the  fit  to  the  3 <>51 
MeV  state o 

A  further  difficulty  arises  in  dealing  with  Z  =  0  capture 
into  levels  of  the  final  nucleus  close  to  the  nucleon  separation  energy 
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FIGURE  4-1-  ANGULAR  DISTRIBUTIONS  FROM  THE  Cl2(d,n)N13 
REACTION  LEADING  TO  THE  GROUND  (J^-  1/2"),  2  367  (JT=l/2f), 
AND  UNSEPARATED  356  (JT  =  5/2*)  *  MEV  STATES  OF  N13  AT  ED  ■  5994  MEV 
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(Ma  60)o  Contrary  to  the  findings  of  Calvert  (Ca  57)  and  Middleton 
(Mi  53)  who  used  9  and  8  MeV  deuterons,  respectively,  the  present  fit 
at  Ep  =  6  MeV  to  the  20367  MeV  level  indicates  £p  =  0  capture  unambig¬ 
uously  o  The  apparent  fit  at  large  angles  is  misleading  and  only 
reflects  the  singularity  of  the  theory  for  that  level 0  At  the  singu¬ 
larity  the  captured  proton  wave  number  is  zero„ 

MacFarlane  and  French  (Ma  60)  suggest  a  simple  Coulomb  cor¬ 
rection  in  the  case  of  (d,n)  transitions  to  such  weakly  bound  final 
states o  This  involves  replacing  the  wave  number  t  of  the  captured 
particle  by 

F2  =  t2  +  A  2M  Ze2  (IV-1) 


where  A  and  M  are  the  target  and  nucleon  masses,  respectively „  This 
correction  was  applied  to  all  the  observed  levels  of  N13„  The  result¬ 
ing  fit  to  the  2 0367  MeV  level  improved  stability  in  that  the  large 
angle  theoretical  distribution  dropped  to  well  below  the  experimental 
distribution,  as  is  customary  for  plane  wave  theory 0  £p  =  0  had  a  more 
reasonable  radius,  r^  =  6  f,  than  £p  =  1  with  7  =>5  f»  The  Coulomb  cor¬ 
rection  did  not  fit  the  unseparated  3 <>56  MeV  level  and  made  no  differ¬ 
ence  to  the  ambiguous  fits  for  the  separated  components  of  the  doublet „ 
The  ground  state  distribution,  unambiguously  Z  =  1,  fitted  better  at 
the  second  maximum  with  inclusion  of  the  Coulomb  effect .  It  should  be 
noted  that  with  increasing  £p  the  stripping  radii  increase,  as  witnessed 
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by  the  third ,  ground*,  and  first  excited  states  having  =  2,  1,  0  and 
rp  =  5  o3  s  and  3  <>6  fermiss  respectively  0 

b0  Distorted  Wave  Analysis 

The  excited  states  of  N^-3  are  unbound  to  proton  emission  and 
therefore  could  not  be  analyzed  by  an  available  distorted  wave  code» 
This  analysis  was  therefore  confined  to  the  ground  state 0  The  only 
distorted  wave  analysis  of  this  state  was  done  by  Hodgson  (Ho  63)  and 
by  Smith  (Sm  63)  with  the  I0U5  to  2 <>95  MeV  data  of  Elwyn,  Kane,  and 
Ofer  (El  59)o  Optical  potential  parameters  were  chosen  consistent 
with  those  found  in  that  work„  The  code  was  one  written  by  Smith 
(Sm  65)  and  previously  adapted  for  the  University  of  Alberta  IBM  70U0 
computer 0  The  optical  potentials  chosen  were  of  the  form 

U(r)  =  Vc(r)  -  V0fws(r,R,a)  -  iUaW  _d  f  (r.R.a)  (IV-1) 

dr 

where  fws(r,R,a)  =  1/(1  =  exp((r  -  rQA1/3)/a)) 

and  Vc(r)  is  the  Coulomb  potential  due  to  a  uniformly  charged  sphere,, 
r  is  the  relative  separation  of  the  reacting  nuclei 9  r0  and  a  are  the 
radius  and  diffuseness  parameters  respectively  for  the  Woods-Saxon 
form  factor;  VQ  and  WQ  are  the  strengths  of  the  refracting  and  surface 
absorbing  parts  of  the  optical  potential 0  The  factor  1*  gives  the 
surface  form  factor  unity  for  its  maximum  value a  The  program  approx¬ 
imates  the  finite  range  n-p  interaction,  the  nonlocality  of  the 
potential,  and  the  deuteron  polarizability . 
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The  parameter  values  are  given  in  Figure  1+-3.  Table  IV-1 

presents  a  list  of  the  optical  potential  parameters  used  for  deuterons 
12 

on  C  o  is  the  reaction  cross  section  and  was  extracted  where  not 
given  and  for  the  present  experiment  from  the  elastic  search  code  of 
Perey  o  The  apparent  violation  of  the  Vr11  ambiguity  in  the  two  sets  of 
parameters  in  the  present  experiment  is  probably  due  to  the  fact  that 
for  deuterons  there  are  discrete  sets  of  VQ  for  a  given  rQ  which  pro¬ 
vide  equally  good  fits  to  the  data  (Pe  63*  Dr  63 ) „ 

The  need  for  a  spin-orbit  term  in  the  incoming  and  outgoing 
channels  is  evident  from  the  too  high  theoretical  distributions  at 
back  angles  (Sh  62).^  Inspection  of  Table  IV-1  reveals  a  real  poten¬ 
tial  approximately  constant  with  bombarding  energy  and  an  imaginary 
term  increasing  with  bombarding  energy 0  This  behaviour  is  consistent 
with  the  investigation  of  Perey  and  Dickens  (Pe  65)0  They  find  that 
the  real  term  is  a  slowly  decreasing  function  of  energy  while  the 
variation  of  the  imaginary  term  with  energy  is  sensitive  to  the  abso¬ 
lute  data  normalization,,  However  the  startling  agreement  with  the 
parameters  of  Smith's  (Sm  63)  2 <,75  MeV  analysis  suggests  that  the 


t  Write  the  scattering  amplitude  as  a  sum  of  a  unit  matrix  and  the 
Pauli  spin  matrices  (for  spin  1/2).  Then  compare  the  two  ratios 
da(7r)/da(0)  =  |a(iT)  |2/|a(0)  |2  with  no  spin-orbit  term  and 
da(Tt)/da(e)  =  |a(7r)  |2/[  |a(e)  |2  +  |  b  ( 0  )  |2sin20]  with  an  £°s_  term. 


where  0  <  0  <  it  0 
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FIGURE  4-3-  ANGULAR  DISTRIBUTION  FROM  C'^d.nJN10  LEADING  TO  THE  GROUND 
STATE  OF  N13-  DISTORTED  WAVE  PREDICTIONS  ARE  SHOWN- 
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Compilation  of  optical  potential  parameters  used  in  12C(d,n)13N  and  12C(d,d)12C  analyses 
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potentials  are  also  sensitive  to  the  code  used0 

Co  J-Dependence  of  the  Angular  Distribution  of  the  3*51  MeV  State 
(J71  =  3/2”)  and  Ground  State  (J^  =  1/2")  of  N13 

Lee  and  Schiffer  (Le  6h )  observed  that  (d»p)  angular  distri¬ 
butions  with  £n  =  1  neutrons  populating  J  =  1/2  states  exhibit  deep 
minima  sometimes  at  about  100°  but  most  often  at  about  130 °0  This 
survey  was  limited  to  mass  number  kO  <  A  <  65  and  energy  7  <_  Ep  <  12 
MeVo  Recently  Davies  et  al0  (Da  66)  found  a  similar  behaviour  still 
perceptible  at  5  MeV  for  the  Si28 (dsn)P29  reaction  leading  to  the 
first  and  second  excited  states  of  P29*  £p  =  20 

The  ground  state  (J77  =  1/2“ )  and  3o51  MeV  state  (J71  =  3/2") 
angular  distribution  are  shown  in  Figure  4-H0  The  apparent  absence 
of  the  Lee-Schiffer  effect  in  this  experiment  should  be  considered 
with  caution  due  to  the  large  error  in  the  angular  distribution  for 
the  second  excited  state  of  N13„  It  is  possible  to  conclude 9  never¬ 
theless*  that  J  dependence  is  probably  absent  at  this  energy  for  this 
nucleus 0  This  is  indicated  by  the  smoothness  of  the  large  angle  dis¬ 
tribution  for  the  ground  state  group 0  The  sparse  data  at  higher 
energies  shows  a  fairly  uniform  large  angle  ground  state  distribution 
as  well  (Ko  63)0 

A  large  angle  comparison  of  the  corresponding  mirror  states  in 
C13  is  not  available 0  The  C12(d*n)C13  ground  state  angular  distri¬ 
bution  is  uniform  between  7  and  11  MeV  (Ev  63)  with  a  dip  at  11  MeV 
and  90  degrees „  Schmidt-Rohr  et  al0  (Sc  63)  however  find  that  the 
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CENTRE  OF  MASS  ANGLE  (DEGREES) 

FIGURE  4-4-  ANGULAR  DISTRIBUTION  OF  THE  GROUND  STATE  VERSUS 
THE  3-51  MEV  STATE  OF  N13-  THE  SMOOTH  CURVES  ARE  DRAWN 
THROUGH  THE  DATA  TO  SHOW  STRUCTURE  IN  THE  ANGULAR  DISTRIBUTIONS 
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ground  state  C13  distribution  is  uniform  at  11 0 8  MeV0 
do  Reduced  Widths 

The  expression  for  the  cross  section  in  a  stripping  reaction 
contains  an  overlap  of  the  initial  target  nucleus  wave  function  with 
the  final  nuclear  wave  function  containing  the  captured  particles  This 
overlap  is  called  the  spectroscopic  factor  S0  If  the  stripped  nucleon 
is  captured  into  a  pure  single  particle  state  so  that  the  passive 
nucleons  are  in  their  initial  state  then  S  =  10  How  well  we  can  deduce 
S  by  comparison  with  experiment  depends  on  how  exactly  we  can  express 
the  remaining  integral  in  the  expression  for  the  cross  sections  This 
contains  the  overlap  of  the  three  wave  functions  describing  the  rela¬ 
tive  motion  of  the  three  active  particles 0 


In  the  DWB  approximation  S  is  related  to  experiment  by 


s  -  °exp 


CIV-2) 


DWB 


In  the  PWB  approximation  with  a  radial  cutoff  a  second  factor,  the 
single  particle  reduced  width  0Q2 ,  appears  as  a  multiplying  factor,, 

The  cross  section  in  this  approximation  is  proportional  to  the  reduced 
width 


( IV-3 ) 


The  general  relationship  (Bu  6l)  is 
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Reduced  widths  were  extracted  in  the  present  experiment  and  are 
tabulated  in  Table  IV-2.  Plane  wave  reduced  widths  were  calculated 
using  formula  (A-17 )  of  French  (Fr  60 ) .  The  value  found  in  this  exper¬ 
iment  for  the  ground  and  3.56  MeV  state  agree  qualitatively  with  those 
found  at  neighboring  energies.  The  singularity  of  the  theory  for  the 
2.367  MeV  state  is  clearly  evident  in  the  large  fluctuations  of  the 
reduced  widths  for  that  state.  The  relatively  low  ground  state  value 
suggests  that  the  target  might  be  thinner  to  the  beam  than  measured. 
Agreement  of  the  plane  wave  spectroscopic  factors  with  those  from  the 
distorted  wave  calculation  would  require  a  single  particle  reduced 
width  different  from  those  found  by  MacFarlane  and  French  (Ma  60 )  for 
C13  by  a  factor  of  five.  This  is  in  agreement  with  the  greater  re¬ 
liability  of  the  DWB  calculation. 

The  DWB  calculation  was  carried  out  with  a  program  using  j-j 
coupling  wave  functions  for  the  final  nucleus,  j-j  coupling  in  fact 
seems  to  be  the  better  choice  for  lp  shell  nuclei  above  mass  10  (La  53), 
although  neither  extreme  coupling  scheme  gives  as  good  agreement  with 
the  level  structure  as  intermediate  coupling  (ibid.,  Co  65)0 
e.  Conclusion 

The  present  experiment  could  not  sufficiently  separate  the 
doublet  in  N1 3  to  present  an  unambiguous  comparison  of  the  ground  and 
second  excited  states.  We  can  however  accept  the  absence  of  the  Lee- 
Schiffer  dip  as  evidence  for  the  lack  of  J  dependence.  This  is 
supported  by  the  lack  of  structure  in  the  ground  state  of  C13  between 
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7  and  11  MeV  (Ev  63,  Sc  63)0  It  would  be  worthwhile  doing  the 
C12(d,p)C13  reaction  in  the  6  MeV  range  as  there  the  doublet  in 
C13  should  be  readily  resolved., 

Bassel  (Ba  62)  points  out  that  the  cross  section  can  be  written 
as  the  sum  of  two  terms 

da  =A(0)  +  (-l)J"jUi/2B(6)  (IV-1) 

dft 

with  A( 0 )  arising  from  the  central  well  distortions  and  B(0)  propor¬ 
tional  to  the  spin-orbit  strengths  of  the  optical  wells ,  In  view  of 
the  conclusion  concerning  the  distorted  wave  fits  in  section  IV  b,  at 
back  angles ,  J  dependence  does  not  seem  to  stem  solely  from  the  spin 
orbit  wells  as  these  are  necessary  to  suppress  the  back  angle  rise  in 


the  theoretical  distribution,, 
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APPENDIX 

a.  General  Derivation  of  the  (d,n)  Cross  Section 

In  a  (d,n)  reaction,  ignoring  center  of  mass  motion,  the  Hamil¬ 
tonian  of  the  colliding  system  can  be  written  in  terms  of  final 
coordinates  (Mo  6 5)  as 

H(£p*£n*C)  =  Ktr^)  +  HF(rp,£)  +  VnI(rn'  ,^)  +  Vnp(  1^'  -  rp|  )  (A-l) 

where  K(rnj=  - (fi2/2MF) Vn2  is  the  kinetic  energy  of  relative  motion, 

Hp  is  the  Hamiltonian  of  the  final  nucleus,  Vnj  is  the  interaction  of 
the  neutron  with  the  nucleons  of  the  target  nucleus  and  Vnp  that  with 
the  captured  proton,,  rp  and  r^'  are  the  coordinates  of  the  captured 
proton  and  of  the  ejected  neutron  relative  to  the  c.  of  m,  of  the 
initial  nucleus,  that  of  the  ejected  neutron  relative  to  the  c0  of 
m.  of  the  final  nucleus  „  Thus  rn'  =  +  [m/(m  +  Mip)]rp.  Mp  is  the 

reduced  mass  for  the  motion  of  the  neutron  relative  to  the  final 
nucleus.  Mp  and  m  are  the  target  and  proton  masses,  respectively. 

The  corresponding  time  independent  Schroedinger  equation  for  the 
system  is 


H<£p. 


In 


—  »rn >£ ) 


(A-2) 


We  next  expand  ¥  in  the  form 
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*(£p.£n.O  =  |^s(£p»^Gs(£n)  (A-3 ) 

where  <fis  is  the  wave  function  for  the  sth  state  of  the  final  nucleus 
of  energy  es ,  so  that 


=  es^s^£p»^  (A-b ) 

and  4>s,  Gs  are  assumed  to  be  complete  in  s0 
Using  relations  (A-l)  and  (A-3),  (A-2)  becomes 


(K(r„) 


+  Hrtrn.C)  +  Vjjjfrn’.C)  +  Vnp( 


±21 


Ep 


Us*. 


(£j,.5)G 


=  E|*s(rp>5)Gs(r11) 

Multiplying  by  and  integrating  over  the  coordinates 

(K<£n)  +  et)Gt(£n)  +  i//d£pd^^t*(In»^)tvnl(£n' 
+  ^pdln'  “  £pl  )]^s(£p,OGs(rn)  =  EGt(rJ1) 


(A-5) 

(£p,4)  yields 
,4) 


or 


-(KfrJ  -  (E  -  =  |Vts  (r„'  )Gs(r„)  (A-6) 


where  E  -  is  the  energy  of  relative  motion  and 


▼*.(£»')  =  JdIt,d5+t*(Ep>5)[vnl(lc'*«)  +  VnpUln'  -  £*,1  ) ]*g (£p.C ) 


(A-7 ) 


Equation  (A-6)  represents  a  set  of  coupled  equations 


xor  une 


wave  function  G^r^)  of  the  relative  motion  in  terms  of  the  set 
4>k(£p,£),  which  describe  the  states  of  the  final  nucleus,  and  the 
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potential  (vnI(rn',5)  +  Vnp(  1^'  -  £p|  )),  Writing  out  in  full  and 
rearranging  slightly  gives 

-(K(rJ  -  (E  -  c0)  +  V00(rn>))G0(rn)  =  0  +  Vqi  W  +  •  .  • 

-(K(£n)  -  (E  -  el>  +  Win-Win)  =  Win’  )Go(In>  +  0  +  •  •  • 

'  '  '  '  (A-8) 

till 

»  »  t  » 

The  set  of  equations  (A-8)  is  general  but  prohibitively  com¬ 
plicated,  Neglect  of  antisymmetrization  with  respect  to  the  exchange 
of  identical  particles  plus  the  following  approximations  may  be 
accounted  for  by  slight  variation  of  the  imaginary  terms  of  the 
optical  potentials. 

In  a  (d#n)  reaction  populating  one  state  k  only,  the  set  (A-8) 
reduces  to 

I 

-(K(£n>  -  <E  -  eo>  +  WiWWln)  "  Win' Win) 

(A-9) 

-(«£n>  -  <E  -  ek>  +  vkk(ln’>Wln>  "  WSn'WStt) 

This  is  the  two  state  coupled  channels  approximation. 

In  the  next  approximation  the  advantage  of  rearrangement  col¬ 
lisions  over  inelastic  scattering  will  become'  evident.  This  is 
due  to  the  fact  that  in  inelastic  scattering,  the  elastic  channel  is 
always  dominant  while  in  rearrangement  there  is  no  source  of  elasti- 
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cally  scattered  outgoing  particles.  I.e.  in  rearrangement  collisions 
^ok^  >  >  ^ko^o  w^ile  the  reverse  inequality  holds  true  in  inelastic 
scattering.  In  this  approximation  we  get  a  pair  of  uncoupled 
equations 


-(KtrJ  -  (E  -  e0)  +  V00(^'  ))o0{r^)  =  Vok(r„’  )%{&) 

,  (A- 10) 

lK(£n)  "  (£!  -  %)  +  )  )G0(rJ1)  ' 

Equation  (A- 10a)  reads 


-(KfrJ  -  (E  -eQ)  +  V00(r„' ))GQ(rn)  (A- 11) 

—  ) E^nl )  +  —  rp I  )]4‘Tr(*'p>£)Gv(rri) 

Since  now  our  expansion  for  V (r^r^jS ) ,  equation  (A-3),  contains 

±.  T- 

only  one  term,  the  k^  1  term,  we  can  equate  <f>k(£p»S  )Gk(£n )  equation 
(A-ll)  to  the  incident  wave 

F0(£d)x(|£0’  -  )ipj(C)  (A-12) 

F^  is  the  wave  function  for  the  deuteron  motion  relative  to  the  initial 
o 

nucleus,  x  is  the  deuteron  wave  function,  is  the  wave  function  of 
the  ground  state  of  the  target  nucleus,  and  r^  =  (r^1  +  jVp)/2,  x  is 
not  to  be  confused  with  the  spin  functions  of  equation  IV-^#  FQ  is  a 
solution  of 


[K(r,j)  -  (E  -  e0 )  +  V^lrj) ]F0(ra)  =  0 


(A-13) 
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3b 


where  eQ  =  0  and 


voo(ld)  =  vdl(ld) 


Fq  is  subject  to  the  boundary  condition 


(A-lU ) 


F0(£a) - *  ♦  MI 


rd  -v  co 


f(0d*^d) 


2? rh2  rd 


(A-ll)  then  becomes 


"(K(£n)  “  (F  -  Go)  +  ^oo^£n ' ) )^o (£n)  “ 

J/VMo',(!i'5,[VnI<In’>«)  +vnp(l£n'  -  rpOH^Oxtlrn'  -  ) 


x  Fo<£d) 


(A-15) 


^©(r^)  can  now  be  solved  for  subject  to  the  boundary  condition 


Go^£n) 


M- 


[F  e 


ik  r 


n  n 


g(  0f*4> ) 


(A-16) 


rn  *  00 


n 


where  k.  2  =  (E  -  e_),  M  is  the  reduced  mass  in  channel  m,  and 

in  — —  m'  *  in  * 


h‘ 


g(0,4>)  is  the  transition  amplitude  for  the  channel  m.  The  solution 
of  (A-15)  can  be  arrived  at  using  the  Green’s  function  method.  First 
the  right  hand  side  of  (A-15)  is  set  equal  to  zero  to  yield  a  homo¬ 
geneous  equation  having  the  same  form  as  the  optical  model  problem. 
Hence  V  can  be  found  from  neutron  scattering  at  the  energy  E,  In 
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inelastic  scattering  we  would  have  to  extract  from  (A-9b )  .  It  is 
impractical  to  observing  elastic  scattering  from  excited  nuclei, 
although  scattering  with  energy  (E  -  e^)  is  sufficient.  At  any  rate 
VQo  is  often  used. 

The  Green's  function  for  (A-15)  can  be  written  in  terms  of 
Go(£n»e )  where 

(K(rJ  -  (E  -  e0 )  +  Vootj^' )  )  =  0  (A-17) 

satisfying 

G0(£n.e)  - *  elknI~n  g(e,*)  (A-l8) 

rn  +  00  2vfi2  rn 

Then 

^(rn)  =  J//'i£ii<irnd5G0(rn ,71-0)4./ (r^.C) 

2-rrh2  rn 

*  tvnl(£n'  +  vnp(  lin'  “  £pl  )^l(5)x(  I £n ’  “  Lp\  )Fo(£d)  (&-19) 

where 

cos  0  =  cos  6  cos  6n  +  sin  0  sin  0n  cos  (4>  -  <J>n)  (A-20) 

At  this  point  it  is  usual  to  neglect  V^r^' ,£ ) .  One  justifica¬ 
tion  for  this  is  to  write 

=  (A-21) 

in  equation  (A-ll)  and  make  use  of  the  orthogonality  property  of  the 
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u)ko  Comparison  of  (A-l6)  with  (A-19)  then  yields 


(A-22) 


and 


da  =  MTMF  kn  g(Qt<p)  2 
dn  (2iTh2)2  kd 


(A-23) 


In  the  foregoing,  spins  have  been  neglected  for  clarity „  In  the  plane 
wave  Born  approximation  the  diagonal  terms  VOG(r^)  and  V00(r^) 9  iee, 
the  optical  plus  Coulomb  potentials,  are  set  equal  to  zero, 
b  0  Programs 

In  the  light  of  the  above  discussion  it  was  felt  necessary  to 
be  able  to  readily  derive  the  potentials  from  elastic  scattering.  The 
recently  operational  magnetic  spectrometer  at  the  University  of  Alberta 
in  fact  provides  a  useful  tool  for  this, 

1 

Two  elastic  search  codes  were  subsequently  converted  for  running 
on  the  University  of  Alberta  IBM  70^0,  These  are  Perey’s  code  and  the 
code  SEEK  by  Melkanoff  (Bu  60 ,  Bu  6l ,  Ma  62,  Me  65*  Me  66,  Ma  66) „  The 
cross  sections  are  found  from  the  scattering  coefficients  C^  which  in 
turn  are  solutions  of 


/  F ^ ' ( p )  +  iGz'(p)  +  C£[F  '(p)  -  iGt'(p)]\ 
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f^  is  the  logarithmic  derivative  determined  from  the  internal  integra- 
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tion0  F^(p)  and  G^(p)  are  the  regular  and  irregular  Coulomb  wave 
functions , 

Perey’s  code  uses  the  Cowell  or  Fox-Goodwin  method  (Fo  ^9 » 
Me  6 5»  Bu  62)  for  the  internal  integration „  For  the  search  the 
function 


x2  =  F(£)  =  l[(q(E)  -  q*)/Afx*]2  (A-25) 

X 

of  the  n  search  parameters  p^ ,  »  .  •  «  Pn  generates  a  true  x2  sur¬ 

face*  f^  and  f^*  are  theoretical  and  experimental  quantities, 
respectively.  The  true  x2”  surface  is  approximated  by  a  quadratic 
X2-  surface  defined  by 


X2  =  F(£)  (A-26) 

where  F(jd)  is  obtained  by  a  second  order  Taylor  expansion  of  F(jd) 
around  the  central  point  jd° 


F(£)  =  F(£°) 
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1  2pi 
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9ZF 


9pi9Pi 


6Pi6p. 


(A-27 ) 


and 


p  =  p  +  6p 


(A-28) 


In  a  neighborhood  of  £°  the  f*  (g)  are  approximated  by  the  first  order 


Taylor  series 


fx(p)  =  fx(p°)  +  1 

X  3Pi 


(A-29) 
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Partial  derivatives  of  F  are  given  by 


3F 

9pi 


=  2Z  (fx<£°)  -  fX*)  3fX 
X  (AfX*)2 


3Pi 


(A-30) 


92F 

2Z  1 

3fx 

3fx 

9pi9pj 

p°  (Af A* )2 

9  Pi 

o  3pj 
p 

(A- 31) 


Substitution  of  eqs »  (A-30)  and  (A-3l)  into  eq.  (A-27 )  yields  the 
quadratic  x2  surface  x2  =  F(p).  Contributions  from  second  partial 
derivatives  of  f^  are  neglected  in  eq.  (A-3l)o  The  first  order  partial 
derivatives  of  f^  are  approximated  by  difference  quotients 

«  fX<Pi°  +  &Pi)  -  fX(Pi°  -  APi} 


Hx 

3pi 


2Api 


requiring  n  +  1  runs  at  p°  and  (p°i,  .  .  PiC„i .  P^°  +  Ap^,  pi°+-.  , 

o  o  o  ,  pn° )  for  all  i o 

Descent  and  termination  proceed  by  the  relation 


3F 

9F 

x  9  92F 

9Pi 

9Pi 

j=l  3Pi3Pj 

P 

P 

6pj  =  0  (i  =  1, 


n) 

(A-33) 


at  the  minimum.  The  solution  of  this  system  of  linear  equation  for  6p 
yields  the  predicted  minimum  point 

p  =  p°  +  6p  (A-28) 


The  further  p  is  from  p°  the  less  reliable  is  the  prediction  of  a  x2 
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minimum,,  Proportionately  smaller  steps  than  6p  are  then  taken  empir¬ 
ically  o  A  parabolic  approximation  with  two  neighboring  values  of  p 
is  used  to  choose  the  next  central  point 0  The  total  number  of  runs 
required  per  central  guess  is  n  +  3.  A  sufficiently  small  improvement 
in  x2  from  one  central  point  to  the  next  central  point  terminates  the 
search  0 

The  code  SEEK  uses  the  modified  Noumerov  (Me  65)  integration 
method  for  the  internal  solution 0  This  is  an  extension  of  the  Fox- 
Goodwin  method  requiring  less  computation  and  no  division „  The  ex¬ 
ploration  phase  of  the  search  is  the  same  as  that  in  Perey's  code 
except  that  the  first  order  partial  derivatives  of  f^  are  calculated 
analytically  using  the  formula  (Me  6U ) 

pm 

_  1  1  1  (4>u)2  dp  u-31*) 

3Pi  ”  E  (on)2  0  3Pi 

oi£  is  a  normalization  constant  and 

V*  =  Vcoul  +  VcR  +  iv=l  +  VSR  +  iVsI 

These  are  the  Coulomb,  complex  central  and  complex  spin-orbit  poten¬ 
tials,  respectively 0  The  derivatives  of  the  cross  sections  in  turn 
are  directly  related  to  3C£/8pio  This  method  of  obtaining  the  quad¬ 
ratic  surface  requires  only  one  run  and  several  numerical  integrations 
as  indicated  by  eq0  (A-3U)  per  central  guess,  This  compares  favour¬ 
ably  with  the  n  +  3  runs  required  by  the  Oak  Ridge  method.  Both 
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codes  require  the  inversion  of  a  p  *  p  matrix  at  each  central  point 
as  indicated  by  equation  (A-33)0 

Perey's  code  was  modified  to  accomodate  a  grid  variation  over 
the  starting  values  of  the  parameters „  Deuteron  polarizability  and  non¬ 
locality  was  included  for  compatability  with  Smith's  DW  (Sm  65,  Di  65) 
code0  Deuteron  polarizability  is  due  to  the  displacement  of  the  charge 
of  the  deuteron  from  its  center  of  mass  and  may  be  accounted  for  by  an 
attractive  dipole  potential 

1  Z2e4ot  (A-36) 

"  2  r4 

external  to  the  nuclear  fieldo  Z  is  the  charge  of  the  target  nucleus „ 
Dickens  and  Perey  (Di  65)  use  a  =  0<,5  (fermis)3  while  Smith  (Sm  65) 
uses  a  =  O0UU5  (fermis )30 

Nonlocality,  which  means  that  the  wave  function  at  r  is  af¬ 
fected  by  its  value  at  r',  can  be  approximated  (Ba  62,  Sm  65)  by 

^NL(r)  =  t1  +  U(r)]”1/2  i|;L(r)  (A-37 ) 

4 

where  U(r)  =  2mV(r)/h2  is  the  real  part  of  the  potential  The  non¬ 
locality  range  3  is  usually  O085  fermis „ 

Table  (A-l)  lists  the  imput  format  to  Perey's  code0  The  poten¬ 
tial  is  of  the  form 

U( r  ,0  )  =  Vc(r,Rc)  -  VQfws(r ,RR,aR)  -cWsfws (r ,Rj ,ax )  (A-38) 

-  iUajWjD  <L_  fws(r >RI>al)  “  £*£.(vso  +  iwso) /_R _  2  i  fws (r  »RR *aR ) 

dr  2mpc  r  dr 
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TABLE  A  -  1 


1 

Card 

Set 

Symbol 

Definition 

Format 

1  ' 

IT 

Function  control 

( 12,18, 

i 

2F10„3,50H 

) 

■ 

norm 

If  ^0,  the  first  NORM  points  are 

0  0  0  0  / 

i 

used  to  search  for  the  best 

1 

normalization. 

j 

WNCUM 

If  #0,  all  experimental  cross 

. 

1 

sections  and  errors  are  multi- 

!  . 

I  ; 

plied,  by  WNCUM  before  calculation 

i 

is  started. 

I  ern 

If  5*0,  the  renormalization  by 

t 

j  ?' 

NORM  is  constrained  from  varying 

I  \ 

too  far  from  unity.  ERN  may  be 

1  i 

taken  as  the  error  in  the 

! 

normalization. 

;  50H 

l  > 

Title 

■ 

2 

AUr 

1 

Read  only  if  17=2,  7  or  8 .  These 

(13F5 ) 

AroR 

are  the  increments  by  which  the 

:  AaR 

(starting)  values  are  changed 

|  AWS 

after  the  completion  of  a  run. 

.  AWq 

AroI 

\ 

A  a  i 

i 

A^so 

i 

i 

AWso 

! 

t 

! 

| 

Ar0c 

! 

! 

A3 

i 

i 

i 

1 

• 

3 

NUGV 

Card  Set  3  must  be  read  in  if 

(11X5) 

NRRGV 

Card  Set  2  is  read  in.  These  are 

: 

i 

nargv 

the  number  of  different  parameter 

\ 

1 

i 

i 

NWSGV 

values  that  each  calculation  is 

i 

nwdgv 

to  be  started  with. 

1 

NRIGV 

NAIGV 

nspgv 

nwpgv 


TABLE  A  -  1  (Continued) 


Card 

Set 

Symbol 

Definition 

Format 

NRCGV 

NBT 

U 

E 

Incident  particle  lab  energy 

(5F10.8) 

AMrp 

Target  Mass 

AM  i 

Incident  particle  mass 

Zt 

Target  charge 

Zl 

Incident  particle  charge 

5 

Ur 

Real  Potential  Well  Depth 

(7F1008) 

roR 

Real  Saxon  well  radius  parameter 

Real  Saxon  well  diffuseness 

parameter 

Ws 

Imaginary  Saxon  well  depth 

WD 

Imaginary  Saxon  derivative  well 

depth 

roI 

Imaginary  Saxon  radius  parameter 

aI 

Imaginary  Saxon  diffuseness 

parameter 

6 

Vso 

Real  Potential  Well  Depth 

(6FIO08) 

wso 

Imaginary  Potential  Well  Depth 

^*os 

Potential  radius  parameter 

a 

Potential  diffuseness  parameter 

b 

3 

Nonlocality  range 

roc 

Coulomb  radius  parameter 

T 

hi 

Integration  step  length  (0„15) 

(F10o8,F5oU, 

a 

Deuteron  polarizability  parameter 

15  ,UF10 08) 

Umax 

Max.,  number  of  partial  waves 

hc 

Coulomb  integration  step  (0ol) 

®min 

Maxo  nOo  of  angles  =  100 „ 

A6 

6  max 

8 

11 

=1  (=  0)  for  spin  (no  spin)  - 

(511) 

orbit  coupling 

12 

=  1  (=  0)  to  print  (not  print) 

(£>»un>JaoO)  X  -  A 


(a.oxne) 
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TABLE  A  -  1  (Continued) 


Card 

Set 

Symbol 

Definition 

Format 

scattering  coefficients,  where 

S^  =  exp(i6^) 

13 

Redundant 

14 

0,  Rj  independent  of  Rr!  1,  Ri=Rr 

15 

0,  aj  independent  of  ap;  1,  aj=aR 

1 6 

1 

Redundant  for  single  calculation , 
For  searchj  0,  do  not  read  in 
o/cjfto  Plot  of  logio(a)  versus  6 

i 

Como 

\  | 

1,  read  a/ap  in  column  31-^0  of 
the  first  NS  cards  of  card  set  10 0 
Plot  of  log10(a/ap)  versus  0  Corn, 

2,  read  c/or  in  column  31-^0  of 
the  first  NS  cards  of  card  set  10 0 

Plot  of  both  logpo(a)  versus  0 

Ccm0  and  logpo (o/ap)  versus  0  c0m9 

j 

9 

Ip  U-2) 

=  0  print  only  final  search 
=  1  print  only  x2 

(12  ,F5.2, 
212,211, 

=  2  print  x2  and  angular 

3I29FlM, 

' 

distribution 

1313) 

AKOFF ( 3-7 ) 

Improvement  in  x2  to  stop  search 

i 

KM( 8-9 ) 

Maxo  no,  of  central  guesses 

1 

allowed  ,  <  100 

+ 

N( 10-11) 

Noo  of  variables  being  searched  on 

; 

* 

Nr(12) 

=  1,  or  searched  on 
=  08  no  search  on  op 

! 

Nt(13) 

=  1,  aT  search  on 
=  0,  no  search  on  a<p 

i 

NS(1U-15) 

NP(l6-lT) 

KR( 18-19) 

no,  of  points  in  aexp(0) 
no,  of  points  in  Pexp(0) 

=  0,  errors  in  a,  P(0)  are  abso- 

! 

i 

i 

lute 

| 

=  1,  A  a,  AP  in  %  of  a(0),  P(0) 

1 

i 

DELM(20-33) 

Min,  value  of  x2  satisfactory 

1 

: 

i 

» 

i 

IPA(I) (3^-36) 

Parameters  being  searched  on 

i 

t]\j  as  o  gives  a  single  calculation  witn  a  pioc 


0  auaiav  (o)qbo I 


lo  4-ic  mst/Xoo  fll  htdi  ts 

-  v  v  (  ■'  [J  ’  ■  - 


?.  .-/tt 


n  ji^.»r  oi  0  = 
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TABLE  A  -  1  (Continued) 


Card 

Set 


Symb  ol 


Definition 


Format 


In  the  current  version  of  the 
program  a  and  3  cannot  be  searched 
upon  o 


10 


0  8a ( 6  )  ,Aa ( 0 )  s  o(0)/oR(0)  if  16=1,2 

not  if  l6  =  0 

0  »P( 0 ) >  AP(0) 

aR>  AaR 
ot»  Aaij 

Ng  +  Np  +  Np  +  Nij  cards  o 


(UF10o8) 


TABLE  A  -  2 


IT 

Function 

Read  Card 
Set  2,3  ? 

0 

Calls  EXITo  This  is  the  last  data  card. 

2 

Single  calculation  with  subsequent  grid-  varia¬ 
tion  of  the  parameters . 

Yes 

3 

Single  calculation  with  no  grid  variation. 

No 

1+ 

Single  search. 

No 

5 

0 

Search  with  same  angles  and  cross  section  as  in 
last  set.  Therefore  omit  card  set  10. 

No 

6 

Search  with  same  data  as  in  last  set.  Read  in 
only  card  set  9. 

No 

T 

Search  with  subsequent  grid  variation  for 
starting  values  of  the  parameters. 

Yes 

8 

As  for  IT  =  T  only  same  angles  and  cross 
sections  as  in  last  set.  Therefore  omit  card 
set  10. 

Yes 

! 

fioi^ooju  i 


,  eTsteuu  iBq  ®rfj  lo  noi* 


i  •  ~  i  a  t  j-  . 


CP  j^a  b*tJ8D  \;Ino 


I  •  iir$  r  .  . 


hi 


where  £-»£  is  for  spin  1/2  particles  only  and  fws(r,R9a)  = 
l/(l+exp((r  -  TqA1 /3  )/a) )  „  Table  (A-2)  contains  the  functions  of 
IT o  All  lengths  and  energies  are  in  fermis  and  MeV,  respectively 0 
c.  Comments  on  the  Use  of  the  Program 

A  six  parameter  search  requires  less  than  10  minutes  of  com¬ 
puter  time.  It  should  be  noted  that  a  search  with  fewer  free  param¬ 
eters  can  often  yield  a  better  fit  as  then  there  are  fewer  local 
minima  in  which  the  search  can  terminate,,  A  need  for  renormalization 
of  the  experimental  data  is  usually  indicated  by  the  preponderance  of 
one  sign  in  the  values  of  chi. 

Perey's  program  is  about  three  times  slower  than  the  code 
Seek.  The  Vrn  ambiguity  can  best  be  treated  by  gridding  over  rG  and 
finding  the  minimum  in  the  U-W-a  space.  In  this  respect  the  grid 
variation  in  the  code  SEEK  is  more  sophisticated  in  that  previous 
best  fit  values  of  the  parameters  are  used  in  subsequent  calculations. 
A  complete  set  of  data  including  polarization  and  reaction  cross 
section  as  well  as  differential  cross  section  generally  removes  such 
amb  iguity . 

The  number  of  partial  waves  available  in  Prey's  code  was  reduced 
from  50  to  HO  1 )  and  the  number  of  experimental  points  from  200  to  100 
in  order  to  calculate  the  Legendre  polynomials  only  once.  This  and 
other  modifications  such  as  the  elimination  of  unnecessary  exponentia¬ 
tion  to  calculate  form  factors  have  speeded  up  the  search  by  about  15$, 


1 )  21.6  MeV  deuterons  on  Mg24  require  a  maximum  l  of  15 . 
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